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Abstract. We present high angular resolution millimeter wavelength continuum and ^^CO(2-l) observations of 
the circumstellar disk surrounding the T Tauri star DG Tauri. We show that the velocity pattern in the inner 
regions of the disk is consistent with Keplerian rotation about a central 0.67 Mq star. The disk rotation is also 
consistent with the toroidal velocity pattern in the initial channel of the optical jet, as inferred from HST spectra 
of the first de-projected 100 AU from the source. Our observations support the tight relationship between disk 
and jet kinematics postulated by the popular magneto-centrifugal models for jet formation and collimation. 
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1. Introduction 

The interplay between accretion and ejection of matter is 
believed to be a crucial element in the formation of stars. 
In particular, stellar jets may contribute substantially to 
the removal of excess angular momentum from the sys- 
tem, thereby allowing the central star to accrete to its fi- 
nal mass ( e.g. E isloffel et al. |2000[ , Konigl & Pudritz pOOC , 
Shu et al. |2000| ). Most of the proposed models invoke the 
simultaneous action of magnetic and centrifugal forces in a 
rotating star/disk system threaded to open magnetic field 
lines. Even if widely accepted, these models have not yet 
been tested observationally on the launching scale (a few 
AU from the star), although this may be possible with 
the coming generation of interferometers. On much larger 
scales, there is some evidence for the requisite relationship 
between the jet and envelope ki nema tics in the HH212 
protostellar system (Davis et al. 2000 ). Hints of rotation 
are seen in the H2 jet knots at 2x10^ to 10* AU from 
the powering source; the sense of rotation is the same as 
that of the flattened envel ope de tected in NH3 VLA ob- 
servations (Wiseman et al. 2001). Although encouraging, 
these measurements probe regions too far from the central 
source to allow detailed comparison between the disk and 
jet kinematics. 

Protostellar disk/jet systems are too embedded to 
probe the jet close to the launching region with current 
techniques which rely on optical and near infrared obser- 



vations. Moreover, the kinematics of disk/envelope sys- 
tems may encompass both rotational and infall motions, 
hampering tests of disk-jet interaction models. Optically 
visible T Tauri stars, which have associated disks but little 
remnant envelopes, are much more suitable candidates for 
such studies. 

The optical jet from the T Tauri star DG Tauri has 
been extensively studied at high resolution in recent years 
and displays properties that are in general agreement with 
magneto-centrifugal models for jet-launching (Dougados 
et al. ^OOOt Bacciotti et al. |2000| ; |2002D . The latter studies 
showed that the flow appears to have an onion-like kine- 
matic structure, with the faster and more collimated flow 
continuously bracketed in a wider and slower one. The flow 
becomes gradually denser and more excited from the edges 
toward the axis. The mass loss rate in the flow is about 
one tenth of the estimated mass accretion rate through the 
disk (Bacciotti et al. 2000). Even more interestingly, for 
the spatially resolved flow component at moderate veloc- 
ity (peaked at — 70kms~^) systematic offsets in the radial 
velocity of the lines have been found in pairs of slits sym- 
metrically located with respect to the jet axis (Bacciotti et 
al. 2002). If these results are interpreted as rotation, then 
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the jet is rotating clockwise (looking toward the source) 
with average toroidal velocities of about 10-15 km s~^, in 
the region probed by the observations (i.e. 10-50 AU from 
both the star and jet axis). All of these properties, includ- 
ing the implied velocities and angular momentum fluxes 
are in the range predicted by the models, assuming a cen- 
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tral star mass of 0.67 Mq (Hartigan et al. 1995). The kine- 
matic properties of the material surrounding the star are 
of crucial importance in further establishing if the models 
apply. 

DG Tauri is known to be surrounded by a circumstel- 
lar disk (Beckwith et al 



199C; Kitamura et al. 1996a 



Dutrey et al. 1996). Previous interferometric observations 



of the molecular component of the system (Sargent & 



Beckwith 1994; Kitamura et al. 1996b, hereafter KKS) 



could not identify a clear signature for rotation around 
the central object. These relatively low-resolution (4-5"), 
^•^00(1-0) observations could not disentangle the kine- 
matics of the circumstellar disk from the outflow and the 
outer envelope velocity fields. In fact, the environment of 
the star on large scales appears to be dominated by out- 
flow motions, possibly due to the interaction between the 
outer regions of the disk and the stellar wind (KKS). In 
contrast, the inner portion of the disk, closer to the jet 
launching region, is expected to display a Keplerian rota- 
tion pattern. We have carried out new, higher-resolution 
millimeter wavelength observations of the ^^C0(2-l) tran- 
sition toward the DG Tauri system with the aim of dis- 
tinguishing the velocity field close to the star and ascer- 
taining if it is consistent with that expected for Keplerian 
rotation and to check that the disk and jet rotate in the 
same sense. 



2. Observations and results 

Millimeter wavelength interferometric observations of the 
DG Tauri system was performed using the Owens Valley 
Radio Observatory (OVRO) mm-array located near Big 
Pine, California, between Oct 1999 and Dec 2001. The 
six 10.4 meter dishes were deployed in configurations that 
provided baselines from 15 to 240 m. Continuum obser- 
vations centered at ~220 and ~108 GHz used an analog 
correlator with a total bandwidth of 2 GHz. The digital 
correlator was configured to observe the ^'^C0(2~l) tran- 
sition with 0.125 MHz resolution over an 8 MHz band 
(0.17 and 11 kms~^, respectively). Frequent observations 
of 0528-1-134 were used to perform phase and gain calibra- 
tion. The passband calibration was obtained by observ- 
ing 3C273, 3C454.3 and/or 3C84. The flux density scale 
was derived by observing Neptune and/or Uranus, and 
the calibration uncertainty is expected to be ~ 20%. All 
calibration and data editing used the MMA software pack- 
age (Scoville et al. 1993| ). Calibrated {u,v) data were then 
loaded into the AIPS and/or GILDAS packages for imag- 
ing, deconvolution and analysis. Continuum maps and line 
cubes were produced using natural weighting of the (u, v) 
data, and smoothed to a spectral resolution of 0.5 kms~^, 
unless specifically noted. The synthesized beam full width 
at half maximum is l'.'7xl'.'4. Continuum subtraction was 
performed on the dirty images before deconvolution using 
channels at the edge of the band. 
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Fig. 1. OVRO 1^00(2-1) integrated intensity map 
(contours) overlaid on the 2.7 mm continuum image 
(greyscale). 



2.1. Continuum maps 

We detect unresolved continuum emission from DG Tauri 
at both 1.3 and 2.7 mm. The peak position is 
the same at both wavelengths, a(2000)=04h27™04!66 
(5(2000)=26°06'16'.'3, in agreement with previous mea- 
surements (e.g. Kitamura et al. 19964 KKS). The to- 
tal flux density is 215 mJy at 222 GHz and 55 mJy at 
108 GHz. Within calibration uncertainties, the 3 mm value 
agrees with earlier interferometer mea surements (KKS; 
Dutrey et al. 1996; Looney et al. 200C) but the 1.3 mm 
value is a facto r of t wo lower than the single dish flux 
(Beckwith et al. 199C ), probably because of spatial filter- 
ing by the interferometer. If we assume optically thin emis- 
sion from dust grains at T~40 K and a dust opacity co- 
efficient ki.=k230GHz X (i^/230 GHz)*^, with k230GHz=0.01 
(Hildebrand 1983 , including a gas to dus t rati o of 100 by 
mass), and f3 ^0.5 (Beckwith & Sargent |l99l| ), our mea- 
surements imply a total mass of '--^0.04 Mq, consistent 
with previous estimates. 



2.2. Line maps 

In Figure [l] we show the ^'^C0(2-l) integrated intensity 
map overlaid on the 2.7 mm continuum image. Our obser- 
vations are sensitive only to the compact features of the 
emission and resolve out most of the extended core emis- 
sion seen in the KKS maps. In spite of the filtering by the 
interferometer, most of the ^^C0(2-l) emission is not con- 
centrated in the inner regions of the system close to the 
position of the optical star. The velocity pattern exhibited 
by this extended gaseous component has been interpreted 
by KKS as due to an expanding disk-like structure, pos- 
sibly the outer edge of the disk that is being dispersed by 
the stellar wind. Our higher angular resolution map is in 
good agreement with their interpretation. In this paper 
however, we will focus on the inner regions of the disk. 
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Fig. 2. ^^C0(2-l) channel maps. Top panels: blue wing; bottom panels: red wing. Each panel is labelled with the 
appropriate velocity {vlsr in kms~^). The cross marks the position of the 1.3 mm continuum peak. Contour levels 
start at 3a and are spaced by la — 60 mJy/beam. The central velocity of the system is assumed to be ~6.0 kms~^ 
(rightmost, isolated panel). 



In Figure ^ we show ^^C0(2-l) channel maps 
(0.5kms~^ resolution) of the central ~12" region centered 
on the continuum peak position (marked by a cross). The 
core of the line (vLSfl=5. 0-7.0 kms~^) is dominated by the 
poorly imaged extended structure discussed above. Note 
that the central velocity channels may also be affected 
by self-absorption due to cold foreground gas (Schuster 
et al. 1993). By contrast, the higher velocity wings, cor- 
responding to the leftmost three blue and red channels, 
display compact emission arising from the inner disk. The 
maps show that the emission peak in the blue channels 
is shifted toward the south-east with respect to the con- 
tinuum peak, while in the red channels it is shifted to 
the north-west. In order to emphasize this velocity gra- 
dient, in Fig. H we show the wing emission integrated 
over the ranges VLSfl=3.5-4.5 kms~^ (blue wing) and 7.5- 
8.5 kms~^ (red wing). These maps were obtained with a 
robust weighting of the (u, v) data and the resulting angu- 
lar resolution is l'.'4xl'.'0. The red and blue wing emission 
peaks on opposite sides of the continuum, which is as- 
sumed to trace the stellar position (see also KKS), and is 
aligned along a line approximately perpendicular to the 
observed direction of the optical jet (p. a. ^226°, marked 
with a thick line in Fig. 0). 
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Fig. 3. '^■^00(2-1) wing maps (contours) overlaid on the 
1.3 mm continuum image (greyscale). Blue and red wing 
are shown as solid and dashed contours, respectively. 
Contour levels start at 3cr and are spaced by la, dotted 
lines show negative contours. The thick solid line indicates 
the direction and extent of the initial segment of the op- 
tical jet studied by Bacciotti et al. (2002). 



3. Discussion 

The prime goal of this study was to investigate the velocity 
pattern in the inner regions of the DG Tauri disk and to 
relate it to the velocity pattern detected at the base of the 
optically visible jet by Bacciotti et al. ( ^002 ). The channel 
and wing maps of Figs. ^ and || indeed show a velocity 
gradient across the inner regions of the circumstellar disk. 
If interpreted as rotation within a disk the axis of which 
coincides with the jet axis, the direction of the gradient is 
consistent with the sense of rotation inferred for the jet. 

Due to the contamination of the poorly imaged exter- 
nal regions of the disk, and self-absorption, it is not pos- 



sible to study the kinematics of the inner disk using the 
line core within 1 kms~^ from the systemic velocity, as- 
sumed to be 5.8 kms~^ (KKS). In particular, there can be 
no detailed comparison of the observed velocity patterns 
with Keplerian rotation models such as those undertaken 
by Koerner et al. (|1993| ), Guilloteau & Dutrey (|1994D and 
Mannings et al. (1997). Nevertheless, in Fig. E position- 



velocity diagrams along directions parallel and perpendic- 
ular to the jet suggest that, along the disk major axis 
(Aj^), higher absolute velocity (with respect to the sys- 
temic velocity) peaks are located closer to the star, while 
the lower absolute velocities peak further away. Moreover, 
blue velocities systematically peak to the south-east and 
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red velocities peak to the north-west of the steUar position. 
This behaviour is quahtatively consistent with Keplerian 
rotation in the inner regions of the disk. A more quantita- 
tive comparison with the expected hne-of-sight velocities 
for a disk surrounding DG Tauri is shown in Fig. |, bot- 
tom panel. The theoretical curves in the figure were com- 
puted for a central star with mass M^,=0.67 Mq (Hartigan 
et al. |1995 ), and an inclination from the line of sight of 
38° (Eisloffel & Mundt |l998|); th ese are the parameters 
adopted by Bacciotti et al. ( 2002 ) to check the rotational 
hypothesis for the jet. The dotted lines include an uncer- 
tainty in these parameters of ±0.25 Mq and ±15°. The 
observed velocity pattern in the inner regions of the disk is 
in excellent agreement with the model predictions. A more 
detailed comparison, including the complete derivation of 
the disk rotation from molecular line observations, will re- 
quire higher angular resolution and more sensitive obser- 
vations of optically thinner transitions, such as C^^0(2-l), 
which are possibly less affected by the external regions of 
the disk/envelope. 

Summarizing our results, we have shown for the first 
time that the disk kinematics in a young T Tauri system 
are in qualitative and quantitative agreement with the ve- 
locity pattern at the base of the jet. In other words, the 
simultaneous and kinematically consistent rotation of disk 
and jet postulated by the popular magneto-centrifugal 
models has been observationally inferred for the first time 
for the region within '^200 AU from the central source. 
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Fig. 4. Top panels: ^"^00(2-1) position-velocity diagrams 
parallel (top) and perpendicular (bottom) to the jet axis. 
The crosses mark the intensity averaged position at each 
velocity, computed only if the emission is above 3cr. 
Errorbars in the top panel show the velocity and spatial 
resolutions. The thin lines mark the core average velocity 
(5.8 kms""'^) and the position of the continuum source. 
Bottom panel: intensity averaged line of sight velocities in 
the first and third quadrant of the lower position-velocity 
diagram, averaged every 0'.'75 (105 AU). The dashed line 
marks the expected line of sight velocity across a Keplerian 
disk with the axis inclined by 38° from the line of sight 
and orbiting about a 0.67 Mq star at a distance of 140 pc 
from the Sun. The dotted lines show the range of variation 
of the Keplerian disk predictions for a maximum uncer- 
tainty of ±15° in the inclination angle and ±0.25 M© in 
the stellar mass. 
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